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ABSTRACT 

During stellar core collapse, which eventually leads to a supernovae explosion, the 
stalled shock is unstable due to the standing accretion shock instability (SASI) . 
This instability induces large-scale non spherical oscillations of the shock, which 
have crucial consequences on the dynamics and the geometry of the explosion. 
While the existence of this instability has been firmly established, its physical 
origin remains somewhat uncertain. Two mechanisms have indeed been proposed 
to explain its linear growth. The first is an advective-acoustic cycle, where the 
instability results from the interplay between advected perturbations (entropy 
and vorticity) and an acoustic wave. The second mechanism is purely acoustic 
and assumes that the shock is able to amplify trapped acoustic waves. Several 
arguments favouring the advective-acoustic cycle have already been proposed, 
however none was entirely conclusive for realistic flow parameters. In this article 
we give two new arguments which unambiguously show that the instability is 
not purely acoustic, and should be attributed to the advective-acoustic cycle. 
First, we extract a radial propagation timescale by comparing the frequencies 
of several unstable harmonics that differ only by their radial structure. The 
extracted time matches the advective-acoustic time but strongly disagrees with 
a purely acoustic interpretation. Second, we present a method to compute purely 
acoustic modes, by artificially removing advected perturbations below the shock. 
All these purely acoustic modes are found to be stable, showing that the advected 
wave is essential to the instability mechanism. 
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1 INTRODUCTION 

Several observational as well as theoretical arguments 
show that most core collapse supernovae are asymmet- 
ric explosions. The observational evide nce include po- 
larisation mea sures jLeonard et al.l l2006). double peaked 



oxygen lines (Mae da et al 



ities of neutron stars (|Lai 



2008), and the large veloc- 
2001). On the theoretical 



side, the failure of realistic sphe rically symmetric calcula - 
tions to produce an explosion l|Liebendorfer et al.ll200ll ) 
suggests that the breaking of the initial spherical sym- 
metry is essential to the explosion mechanism. Indeed 
most proposed mechanisms are fundamentally asymmet- 
ric: th e modern version of the neutrino-driven mech- 
anism (|Marek fc Jankal 1200 9*). the acoustic mechanism 
( Bur rows et al ■1120061). and t he magneto-centrifugal mech- 
anism (| Wheeler et al.ll2002l ). 



The spherical symmetry can be broken by several 
phenomena. In the magneto-centrifugal mechanism, the 
asymmetry is due to a fast rotation in conjunction with 
a strong magnetic field. In the case of moderate ro- 
tation, two non radial instabilities can develop during 
the stalled shock phase: the neutrino-driven convection, 
and the stationary accretion shock instability usually 
called SAS I which induces non spherical oscillations of 
the shock |Blondin et al.ll2003l ). Of these two instabili- 
ties, S ASI is the most efficie nt in creating a global asym- 



metrv (|Foglizzo et al, 20061) and helping t he explosion 
(|Marek fc Jankal 12009 : iBurrows et al.ll2006l ). In addition 



to its crucial role in the explosion mechanism, SASI has a 
number of potentially important consequen ces: it affects 
the kick and the spin of the neutro n star l|Scheck et al.l 
120061 : iBlondin fc Mezzacappall2007l ). causes an emission 
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of gravitational waves (|Kotake et al J 120071 1. and creates 
a tim e dependence of the neutrino signal (jMarek et al.l 
l2009h . 

But what is SASI ? Which physical phenomenon 
drives these shock oscillations ? Surprisingly, the an- 
swer to this basic question remains somewhat uncer- 
tain. Indeed, two mechanis ms have been propo sed - 
the advective-acoustic cycle (Foglizzo ct al. 2007) and 
a pu rely acoustic mechanism (jBlondin fc Mezzacappal 
120061 ) - and distinguishing them remains a controversial 
issue. The aim of this article is to clarify this question. 

In Section[2] we review previous works by explaining 
the two instability mechanisms and the arguments pro- 
posed to distinguish them. On the basis of these works, 
the advective-acoustic cycle is generally favoured but the 
purely acoustic mechanism cannot be ruled out for re- 
alistic flow parameters. In Section [3] we propose a new 
method to extract a radial propagation time from the fre- 
quency of the unstable modes. For this purpose we com- 
pare the frequencies of several harmonics that differ only 
by their radial structure. As the radial advective-acoustic 
and purely acoustic times are significantly different, this 
allows us to distinguish between the two cycles and dis- 
cards the purely acoustic mechanism. In Section^ we de- 
scribe a method to compute purely acoustic modes, which 
are shown to be stable. In Section [5] we give a detailed 
analysis of the frequency spectra of advective-acoustic 
and purely acoustic cycles, which confirms the validity 
of the method used in Section [3] Finally in Section [6] 
we summarise our work and conclude that the instability 
mechanism behind SASI is the advective-acoustic cycle. 



2 THE TWO PROPOSED MECHANISMS: 
SUMMARY OF PREVIOUS WORKS 

2.1 Different models of SASI 

SASI has been observed to grow in many different 
models of stellar collapse, which have very diverse de- 
grees of complexity. These incl ude state of the art 
realistic numerical simulations ( Ma rek fc Jankal 120091 ; 
I Burrows et al.l 120061 ') , as well as models where succes- 
sive simplificati ons have been mad e: approximate neu- 
trino transport (Schcck ct al. 2008 ), neglect of the stel- 
lar structure (jOhnishi et al.l 120061 ; lYamasaki fc Yam ada 
120071 ; llwakami et al. 20081 ). simple equation of state with 
the effect of neutrinos pa r ameterised by cooling/heating 
funct i ons jBlondin et al] 120031; iBlondin fc Mezzacappal 
20061; IFoglizzo et alj 120071; lYamasaki fc Foglizzd |200S| ; 
Fernandez fc Thompson 2009b. a), adiabatic approxima- 



tion (|Blondin fc Mezzacappall2007l ). and finally the most 
simplified of a ll these models is the planar toy model 
(|Foglizzd 120091 ; ISato et al.ll2009l ). Very often, the sim- 
pler the model, the deeper the physical understanding. 
For example, models neglecting the heating have clearly 
show n that SASI is distin ct from neutrino-driven convec- 
tion l|Blondin et al.l 12003). a fact missed by most realis- 
tic models. The neglect of the stellar structure allows 
to set up an initial steady state, which can be stud- 
ied with a perturbative a nalysis (|Foglizzo et al. 1 120071 ; 
lYamasaki fc Yamadall2007l ). Finally, the simplest of these 



models, the planar toy model, allows for an analytical 
treatment where the instability mechanism can be clearly 
demonstrated. 

It is widely accepted that the same instabil- 
ity is at work in all these models. Therefore, al- 
though we base our i nvesti gation on the model of 
IBlondin fc Mezzacappal (j^OOS) , we expect our conclu- 
sions to hold more generally. We choose this model 
because it is simple enough to perform a perturba- 
tive analysis with clean boundary c onditions. This ide- 
alised flow was first described by 



1992) 



Houck fc Chevalier 



and later studied by IBlondin fc Mezzacappa 
20061 ) to describe SASI. It has then become a clas- 
sical model of SASI and has been studied also by 



Foglizzo et all J2007h; iFernandez fc Thompson! l|2009bl ); 
Fernandez! j2010l );l Guilet et al. (|2010l ). In this flow, the 



fluid is modelled as a non self-gravitating perfect gas with 
an adiabatic index of 4/3. The effects of neutrinos are 
simply modelled by a cooling function, while neutrino 
heating is neglected. The proto-neutron star is treated as 
a hard surface on which the cooling fluid settles down. 
The shock is stationary at a radius which can be freely 
chosen by adjusting the normalisation of the cooling func- 
tion. This model and its linear analysis will be extensively 
us ed in the remainder o f the paper. We refer the reader 
to IFoglizzo et all (|2007l ) for a description of the linear 
eigenmodes calculation, and to Appendix|X]for a slightly 
different (but equivalent) formulation of the equations 
governing it, which has been used in this paper. 



2.2 The two mechanisms 

The two mechanisms proposed to explain the growth 
of SASI are schematically illustrated by Figure [1] The 
advective-acoustic cycle (left panel) is a cycle between 
two waves : an advected wave composed of entropy 
and vorticity and an acoustic wave propagating outward 
(IFoglizzo fc Tagger! l2000l; iFoglizzol [20021 ; IBlondin et all 



120031 ; Foglizzo et al.1 20071 ). Note that the acoustic wave 



also propagates in the transverse direction. These waves 
are coupled via two coupling processes taking place at the 
shock and at a smaller radius in the deceleration region. 
The acoustic wave reaching the shock makes it oscillate 
and creates the entropy-vorticity wave. This wave is ad- 
vected towards the proto-neutron star (hereafter called 
PNS) and its deceleration creates an acoustic feedback 
thus closing the loop. The cycle is unstable if the ampli- 
tude of the wave at the end of a cycle is larger than it 
was at the beginning. 

The purely acoustic m e chani sm has been described 
in IBlondin fc Mezzacappal (|2006l ) as an acoustic wave 
trapped in the postshock cavity. This acoustic mode is 
assumed to be amplied through its interaction with the 
shock. Even when the acoustic propagation is mostly az- 
imuthal some radial propagation is also present between 
the shock wave and th e inne r turning point (Fig. 6 of 
IBlondin fc Mezzacappal (2006)). Such an acoustic path 
can be interpreted as a purely acoustic cycle between the 
shock and the turning point (Figure [TJ right panel). We 
stress that this description is quite general and does not 
assume a purely radial propagation. The transverse prop- 
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Figure 1. Schematic representation of the linear instability mechanism proposed to explain SASI (adapted from [Foglizzo et al.l 
| |2009| )). The advective-acoustic cycle (left panel) comprises two coupling processes between an advected wave (represented here 
by the circular red arrows) and an upward propagating acoustic wave (represented by wavy arrows) . These coupling processes take 
place at the shock and in a region of strong deceleration close to the surface of the proto-neutron star. Note that while advection 
is purely radial, acoustic propagation can have a significant transverse component. In the purely acoustic mechanism a trapped 
acoustic wave is amplified through its interaction with the shock. Formally, it can be represented by two coupling processes between 
two acoustic waves, propagating up and down (right panel, note the concurrent transverse propagation). The advective-acoustic 
and purely acoustic cycles have been represented respectively unstable and stable, in anticipation of the conclusion of this article. 



agation is not restricted or disregarded neither here nor 
anywhere in this paper. 

Any of these two cycles could in principle be un- 
stable depending on the efficiency of the different cou- 
pling processes. For example, the purely acoustic cy- 
cle might be unstable if the acoustic wave is ampli- 
fied _wheji_Jt_^s_^flected a t the shock, as assumed 
by iBlondin fc Mezzacappal |2006). Alternatively, the 
advective-acoustic cycle can be unstable if gradients be- 
low the shock allow for an efficient acoustic feedback, 
and if the shock is stro ng enough to efficiently create an 
entropy- vo r ticity wave (|Foglizzo fc Taggeill2000l ; iFoglizzol 
l200ll . 120021 , |2009h . By stu dying a s i mplifi ed model of the 
advective-acoustic cycle, IFoglizzol (12009 ) could explain 
salient features of SASI, in particular that it is dominated 
by low frequency and large scale non-radial modes. On 
the other hand, the purely acoustic cycle lacks a reference 
model where the properties of such acoustic instability 
could be unambiguously determined. As a consequence, 
it has not so far explained the basic features of SASI. 



2.3 Timescales 

A large part of this article deals with timescales: in par- 
ticular how they determine the oscillation frequency, and 
how this can be used to discriminate between the two 
mechanisms. It is therefore useful to recall the character- 
istic timescales of the two cycles. The duration of a radial 
advective acoustic cycle can be written as: 

1 dr 



precise location of this coupling. Similarly, the duration 
of a radial purely acoustic cycle is: 



• MM' 



(1) 



where r CO up is the effective c oupling radius where the 
acoustic feed b ack is producecFl lFoglizzo et al.l (|2007h and 
IScheck et alj (|2008l ) have suggested that the coupling ra- 
dius is close to the location where the deceleration is 
strongest, however there remains some uncertainty in the 



1 Note that this is a simplified description of the instability: 
in reality the feedback is distributed radially rather than pro- 
duced at a single radius. 



dr 



M 2 



(2) 



where r rc f is the radius where the acoustic wave propa- 
gating down is reflected (i.e. its turning point). Finally, 
another acoustic timescale which plays a role potentially 
for both of the cycles is the azimuthal acoustic time: 

27rr 



(3) 



The two radial timescales suffer some uncertainty due to 
the lack of knowledge about r coup and r re f , both of which 
can in principle depend non trivially on the mode con- 
sidered. In the model studied in this paper and briefly 
described in Section \2. II it is however possible to obtain 
a useful upper limit by supposing that these two radii co- 
incide with the PNS surfac^E In the following, we choose 
this prescription, keeping in mind that a shorter time 
(hence a higher frequency) can be obtained by changing 
r COU p and r re f. The azimuthal acoustic time also suffers 
an uncertainty due to the radius at which it should be es- 
timated. Here again a useful upper limit can be obtained 
by choosing the shock radius r sh . 



2.4 WKB des cription of the t wo cycles 
(following iFoglizzo et al.l (120071 )) 

In order to quantitatively describe the two cycles, one 
needs to separate perturbations below the shock into the 
three types of waves: two acoustic waves (propagating 
up and down) and the advec ted wave. Below the shock 
the flow is close to adiabatic (|Foglizzo et al l (|2007l ). and 



2 With the cooling function considered here (a = 3/2, f} = 
5/2), the advection time until the PNS surface is finite. How- 
ever with a less violent cooling law (a = 1, j3 = 6) this time 
is infinite, so that no useful upper limit can be obtained. It 
is then more relevant to estimate the coupling radius to b e 
where the velocity gradient is strongest llFoglizzo et aljfeoOTl) . 
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see Section [4]), thus the entropy and vorticity are sim- 
ply advected. The main difficulty comes from the acous- 
tic perturbations, which can be separated between waves 
propagating up and down only in the WKB approxima- 
tion. The WKB approximation is valid if the wavelengths 
of the acoustic and advected waves are small compared 
to the length scale of the gradients. It can therefore be 
applied only to the modes with a high enough frequency, 
i.e. satisfying: 



by: 



log IQI 



w — > 1. 

Csh 



(4) 



iFoglizzo et al] (2007) performed such a WKB a nalys is of 
the flow considered bv lBlondin fc Mezzacappal l|200rJ ). In 
the following, we summarise the principle of this analysis, 
which could in principle be applied more generally to 
other flows. 

Constants describing the coupling efficiency can be 
defined in the following way. When an acoustic wave 
propagating up with an amplitude 8f~ hits the shock, 
it creates an advected wave with an amplitude <5/ adv as 
well as an acoustic wave propagating down with an am- 
plitude 8f+ 0. The ratio of the amplitude of the waves 
created to that of the incoming wave defines the coupling 
constants: Q sh = 5/ adv /<5/ _ and R sh = 5f + /Sf~. These 
can be computed analytically using the boundary condi- 
tions at the shock and by separating the perturbations 
into the three types of waves owing to the WKB method 



into the three types ot waves owing t 
l|Foglizzo et al.ll2005l : IFogl izzo 2009). 



When an entropy- vorticity wave is advected towards 
the PNS, it creates an acoustic wave propagating up 
with an efficiency characterised by the ratio of the ampli- 
tudes (measured at the shock): Qv = 5/~/<5/ adv . Simi- 
larly the reflection in the gradients of the acoustic wave 
propagating down is characterised by: i?v = Sf~ /Sf + . 
These two constants can be computed numerically by in- 
tegrating the differential system describing the pertur- 
bations, imposing modified boundary conditions at the 
shock (Sf + = to obtain Q v , <5/ adv = to obtain R v ). 

Finally, using these four coupling constants, one can 
define the global constants that describe the efficiency 
of the two cycles (the amplitude ratio between the be- 
ginning and the end of the cycles): Q = Q B hQv and 
R = _R s h_Rv- The cycle constants are complex numbers 
that depend on the complex frequency u. Since the two 
cycles take place simultane ously, an eigenmode mus t sat- 
isfy the following relation (jFoglizzo fc Tagger 20o3): 



' + R= 1. 



(5) 



If these constants are calculated with a real fre- 
quency then the instability criterion for a single cycle can 
be written: \Q\ > 1 or \R\ > 1 (depending on the cycle 
considered) , which is equivalent to the criterion stated in 
Section [2.21 The growth rate can then be approximated 



3 the definition of the variable Sf is recalled in Appendix 1X1 
Choosing another variable to measure the amplitude would 
change the coupling constants at the shock and in the gradi- 
ents, but not the global cycle constant determining the stabil- 
ity of the cycles. 



(6) 



where r is the duration of the cycle considered (here the 
advective-acoustic one). When a second cycle is present, 
it can interfere constructively or destructively with the 
first cycle depending on their relative phase. This causes 
oscillations of the growth rate, which are i ndeed observed 
in the eigenspectrum of SASI (Figure 7 of Fogl izzo et 
(|2007h ). These oscillations contain information about the 
efficiency and the t imescale of the two cycles at work. 
IFoglizzo etHI l|2007h could thus extract directly from the 
eigenspectrum the values of the constants Q and R (note 
that this is independent from the WKB analysis, and 
only assumes the presence of two cycles). The value of 
the cycle constants obtained from the WKB method and 
from the analysis of the eigenspectrum are in very good 
agreement (their Figures 8 and 9). They show that the 
advective-acoustic can be unstable with a cycle constant 
reaching Q ~ 6, while the purely acoustic cycle is always 
stable with typically R ~ 0.5. 



2.5 Which cycle drives SASI? 

The WKB analysis of IFoglizzo et al.l <|2007h summarised 
in the last subsection has proven that the advective- 
acoustic cycle is unstable and can explain SASI, while the 
purely acoustic cycle is stable. However this result is valid 
only for high frequency modes, satisfying Equation [4] 
The first few harmonics do not satisfy this inequality as 
o;r s h/c s h ~ 1 for the fundamental mode. One may con- 
sider that the WKB method can be safely applied for the 
tenth harmonics and higher (satisfying wr s i,/c s h > 10). 
Unfortunately these high frequency modes are unstable 
only for large shock radii r s h/r» > 5 — 10 (depending on 
the cooling function chosen). Thus strictly speaking in 
the more realistic range of shock radius 2 < r s h/r» < 5, 
the WKB analysis cannot concl ude on the instability 
mechanism. IFoglizzo et all <|2007h have however argued 
by a continuity argument that it is not necessary to in- 
voke a different instability mechanism for the low fre- 
quency modes. Indeed, the eigenspectrum vary smoothly 
both with the shock radius and mode frequency, and the 
eigenfunction of the low frequency modes resembles that 
of the high frequency modes. It may thus seem natural 
to assume that all SASI modes are due to the advective- 
acoustic cycle. Another argument in this direction is that 
the WKB approximation usually gives quite good results 
even when the small parameter (here c s h/(r a hu;)) is of or- 
der unity. This is however not a formal proof, and it can- 
not be excluded that the low frequency modes at mod- 
erate shock radius originate from a different instability 

mec hanism. 

iLamind l|2007h has studied analytically the stability 
of a spherical accretion shock by establishing an approxi- 
mate dispersion r elatio n, using a method inspired from 
IVishniac fc Rvul (|l9891 ). This equation contains terms 
due to advection, which lEami ng (2007) proposes to re- 
move in order to assess the importance of the advective- 
acoustic cycle in the instability mechanism. The conclu- 
sion he initially reached by this mean depended on the 
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radius of the shock. At large shock radii, the omission 
of the advective term stabilised the flow, thus favour- 
ing the advective-acoustic cycle. On the contrary, at 
small shock radii advection affected very little the growth 
rate, favouring a purely acoustic interpretation of the in- 
stability. However, this ana lysis has been criticised by 
lYam asaki fc Foglizzol (|2008l ) because the dispersion re- 
lation was established with the use of dubious approx- 
imations. FWthermore an erratum published later on 
(|Lamindl2008h corrects a few mistakes and changes the 
conclusion: the advective-acoustic cycle is then favoured 
for all shock radii. 

A number of authors have tried to distinguish the 
two mechanisms by using the oscillation frequency ob- 
served in the simul ations. These studies have led to di- 
verging conclusions: lOhnishi et alJ (|2006h concluded that 
the advecti ve-acoustic cycle is responsible for the insta- 
bility, while iBlondin fc Mezzacappal (I2006T I concluded in 
the favour of the purely acoustic mechanism. This confu- 
sion probably comes from the fact that both of the mech- 
anisms can explain the oscillati on frequency if one m akes 
the right assumptions. Indeed, lOhnishi et alJ (|2006l ) ex- 
cluded the purely acoustic mechanism because the os- 
cillation per iod is much longer than the r adial acoustic 
time, while IBlondin fc Mezzacappal l)200rJ ) have shown 
that the oscillation period could be explained if one con- 
siders a non radial acoustic p ath. On the other hand, 
IBlondin fc Mezzacappal (|2006) excluded the advective- 
acoustic cycle without considering the possibility that 
the acoustic feedback be created at a larger radius than 
the PNS surface. The determination of the frequency of 
each cycle thus suffers some uncertainty due to our lack 
of knowledge of the coupling radius (for the advective- 
acoustic cycle), and of the non-radial pat h to be con- 
sidere d (for a purely acoustic cycle). Later. [Scheck et al.l 
(|2008T ) have shown that the oscillation period observed 
in their simulations was compatible with both mecha- 
nisms if the advective-acoustic coupling radius lies where 
the deceleration is maximum, or if an ad hoc non radial 
aco ustic path is chosen (their Figure 16). 

iFernandez fc Thompson! (|2009b) have proposed an 
argument that can be applied to the regime in which the 
WKB method is inapplicable. They followed the varia- 
tion of a SASI mode growth rate as the shock radius is 
varied, which changes the ratio of the azimuthal acoustic 
time to the advection time from the shock to the PNS. 
They observed that the mode has a maximum growth 
rate when the two times coincide, and interpreted this 
as an indication that the advective-acoustic cycle is at 
work (because both advection and acoustic times play a 
role in the instability). While this observation is indeed 
suggestive that the advection plays a role in determining 
the growth rate, it does not rule out the possibility of the 
acoustic cycle playing a dominant role with only a minor 
help from the advection. 

Some more indirect information about the instabil- 
ity mechanism can be found by considering the effect 
of a mag netic field on SASI. lEndeve et all (|2010h have 
performed numerical simulations of SASI in the presence 
of a radial magnetic field (a split monopole). They ob- 
tained the surprising result that a magnetic field, which 
magnetic pressure is negligible compared the thermal one 
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Figure 2. Comparison of the oscillation period of the fun- 
damental I = 1 SASI mode with the advective-acoustic and 
purely acoustic times. This calculation is made in the classi- 
cal model of SASI briefly described in Section 12.11 with the 
parametters a = 3/2, /3 = 5/2, M\ = co (as in the rest of 
the paper). All timescales are in units of r^/v^. The oscil- 
lation period (t OS c, red line) is comparable both to the radial 
advective-acoustic time (t aac , full black line, Equation JT)) 
and to the azimuthal acoustic time at the shock (t ac ,*, short 
dashed line, Equation 101). By contrast, the radial acoustic 
time is much shorter than these timescales (t aC r, long dashed 
line, Equation J2j). For all timescales an upper limit is shown 
as discussed in Section 12.31 (r C oup = »"*, r rc f = r„, and r ac is 
estimated at the shock radius). 

(fmag/fth < 10~ 3 ), is capable of entirely stabilising SASI 
(in their model 2DB14Am). This seems to contradict a 
purely acoustic mechanism because the acoustic waves 
should be almost unaffected by such a weak magnetic 
field. By contrast, this result can be naturally interpreted 
in the framework of the advective-acoustic cycle. Indeed 
iGuilet fc Foglizzol (|201Cj ) have shown that the advective- 
acoustic cycle is significantly affected by a magnetic field 
if the Alfven speed is comparable to the advection ve- 
locity, which can happen even if the magnetic pressure 
is negligibly small. In their model 2DB14Am, the Alfven 
speed is indeed comparable to the advection velocity in 
the vicinity of the PNS. 

To summarise this subsection, the advective-acoustic 
mechanism is favoured by most arguments. However none 
of these arguments can rule out the purely acoustic mech- 
anism for realistic flow parameters. The goal of this pa- 
per is to find a method to definitely distinguish the two 
mechanisms for realistic flow parameters. 

2.6 This article 

As discussed in the last subsection, previous studies were 
not successful in determining the instability mechanism 
from SASI oscillation frequency, because the fundamen- 
tal mode frequency can be explained by both mecha- 
nisms. This is due to the fact, that the (radial) advective- 
acoustic time and the azimuthal acoustic time are both 
comparable to the oscillation period, as is illustrated by 
Figure [2] In Section [3] we present a method to avoid this 
difficulty by extracting a radial time from the frequencies 
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of several unstable radial harmonics. The large difference 
between acoustic and advective-acoustic radial times al- 
lows us to distinguish clearly the two mechanisms (the 
radial acoustic time is shorter by a factor 3 to 10, as 
shown in Figure [2}. 

Another reason for the remaining ambiguity on the 
instability mechanism is that, for realistic flow param- 
eters, the unstable modes cannot be described by the 
WKB method. In Section [4] we present a new method to 
compute purely acoustic modes, which partially avoids 
this difficulty. This method makes use of the fact that the 
wavelength of advected perturbations is shorter than the 
acoustic wavelength. As a consequence, the identification 
of advected perturbations through a WKB approxima- 
tion is valid at lower frequency. 

Finally, in Section [5] we develop an analytical de- 
scription of the frequency spectrum of the two cycles, in 
order to confirm the validity of the argument presented 
in Section[3l This analytical treatment is checked against 
the acoustic modes computed in Section 2] 

In the following we apply our two arguments to the 
classical model of SASI briefly described in Section 12.11 
The reader is referred tolBlondin fc Mezzacappal (|200fJ ) 
an d lFoglizzo etafl (|2007h for a more detailed description. 
We focus on a shock radius r s h = 2.5r», which is typical 
of supernovae and for which the unstable modes cannot 
be described by the WKB method. We also choose the 
following parameters: a — 3/2, /3 = 5/2, Mi = oo. 



3 EXTRACTION OF A RADIAL 

PROPAGATION TIMESCALE FROM THE 
EIGENSPECTRUM 

Previous studies that tried to determine the in- 
stability mechanism us i ng th e oscillation fre q uency 



Rectangular box 



stability mechanism us i ng tri e oscillation fre q uency 
(IBlondin fc Mezzacappal |2006| ; lOhnishi et all 12009 ; 
IScheck et al. 20081 ) have based their discussion entirely 
on the mode that grew fastest in their simulation. This 
most unstable mode is usually the lowest frequency 
I = 1 mode, and its frequency is consistent with a radial 
advective-acoustic time as well as with a non-radial 
acoustic time. The specificity of our study is that 
we consider other (less) unstable modes, which can 
provide us with additional constraints on the instability 
mechanism. Th ese modes are accura tely computed by a 
linear analysis (|Foglizzo et all 120071 ). but may be hard 
to distinguish clearly in numerical simulations (although 
they should be present) because they grow more slowly 
and thus have a smaller amplitude than the most un- 
stable mode. We focus in particular on higher frequency 
harmonics having the same transverse structure as the 
most unstable 1 = 1 mode. For most astrophysically 
relevant values of the shock radius, several of these 
modes are unstable. For example, for our fiducial value 
r s h = 2.5r«, three I = 1 modes are unstable (Figure [10}. 
We emphasise that these modes have a higher frequency 
than the fundamental mode, but not high enough to 
allow for a WKB analysis. 

What timescale determines the frequency change of 
the higher harmonics as compared to the fundamen- 
tal mode ? We argue that this timescale cannot cor- 



SASI cavity 



Shock 




Figure 4. Schematic illustration of the analogy between 
acoustic modes in a rectangular box and in stellar core col- 
lapse. 



respond to an azimuthal propagation because all these 
modes share the same transverse structure (determined 
by the spherical harmonics index I). As illustrated by 
Figure [3] these modes differ only by their radial struc- 
ture suggesting that the frequency change should instead 
be controlled by a radial time, either advective-acoustic 
or purely acoustic. As a consequence, we suggest that a 
radial propagation timescale can be extracted from the 
frequency difference between two successive modes that 
share the same transverse structure (i.e. same spherical 
harmonics). We define such a time in the following way: 



trt- 



2tt 



u) r (n + 1) — ui r (n) 



(7) 



where uv(n) and uj r (n + 1) are the frequencies of the n-th 
and n+l-th mode of SASI for a given spherical harmonics 
index I (ranking the modes from the lowest to the highest 
frequency). From the physical argument stated above, we 
expect this time to be of the order of the radial acoustic 
time if the instability mechanism were purely acoustic, 
or of the order of the radial advective-acoustic time if 
SASI were due to the advective-acoustic cycle. Given the 
qualitative nature of the argument, we may not however 
expect a precise agreement. Indeed the analysis of the fre- 
quency spectrum detailed in Section [5] confirms that this 
is only approximately true. However thanks to the large 
difference between radial times (Figure [2}, this approxi- 
mate estimate is enough to be conclusive. We note that 
this method bears some similarity with the extraction of 
the radial structure of the S un using the "large freque ncy 
spacing" in helioseismology l|Deubner fc Gougbll 19841 ). 

Some light can be shed on the above argument by 
drawing an analogy with the very classical case of acous- 
tic modes in a rectangular box (illustrated by Figure [4}. 
Let us assume that the box is filled with a uniform gas at 
rest and delimited by reflecting rigid walls. The frequency 
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x x 



Figure 3. Normalised pressure perturbations 8P/P in the equatorial plane associated to different SASI modes. I = 1, m = 1 modes 
are depicted in the left column, while I = 2, m = 2 are in the right column. The three panels in each column correspond to the 
three lowest frequency modes for the given spherical harmonics index (low to high frequency represented from top to bottom) . This 
illustrates how modes with the same spherical harmonics index differ by their radial structure while they share their transverse 
structure: higher frequency harmonics have a more complex radial structure. 
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eigenspectrum of acoustic modes is then: 



= 2tt 
= 2tt- 



2 nl/2 



1 + 



1/2 



(8) 
(9) 



where Tacz = 2L x /c and r aC z = 2L z /c are the back and 
forth acoustic propagation times in the x and z direc- 
tions, while n x and n z are the number of half wave- 
lengths in the width of the box (in x and z direction). 
To draw the analogy, we consider a box elongated in 
the ^-direction, with z representing the radial direction 
and x representing the transverse direction (8 or (f>) . The 
correspondences are thus L x — > nr, L z — > r s h — r* and 
n x — > I. n z distinguishes the different modes that share 
the same spherical harmonics index. n z = corresponds 
to the lowest frequency fundamental mode, with a fre- 
quency uj = 2irn x /T a , cx that is dictated by the transverse 
acoustic time. Higher integer values of n z correspond to 
the higher frequency modes, which are the focus of this 
section. As n z is increased the frequency increases due to 
the first term in Equation JSJ, which contains the acous- 
tic time Tacz in the ^-direction (analogous to the radial 
acoustic time). This confirms the qualitative argument 
given in the last paragraph. As the frequency is not sim- 
ply proportional to 2nn z /T !!lcz , the frequency difference 
between two successive modes is not exactly the "radial" 
frequency 2%/T acz , suggesting that the Equation Q is 
only an approximate estimate even in this simplest case. 
However the extracted timescale does tend toward the 
radial time in the limit where n x /t£ cx <C n 2 z /r acz . As the 
azimuthal acoustic time is significantly larger the radial 
time, we can expect the method of extraction of a radial 
timescale to give reasonable results. 

The radial timescales extracted by using the mode 
duplets n = 1 — 2, n — 2 — 3, and n — 3 — 4 (all of them 
with 1 = 1) are compared to the advective-acoustic and 
purely acoustic radial times in Figure \S\ The timescales 
extracted from mode duplets are represented only on the 
radius interval in which both of the modes are unstable. 
The three extracted timescales are close to each other, 
which gives some confidence on their physical relevance. 
Furthermore, they roughly match the advective-acoustic 
time: within ~ 15% for moderate shock radius in the 
range r s h = 2 — 5r» (the focus of this article), and within 
~ 50% for larger shock radii. By contrast, the extracted 
timescales are significantly longer than the radial acoustic 
time: the ratio varies between 3.5 for small shock radii 
and 6 for larger shock radii. This result clearly excludes 
the purely acoustic mechanism and argues in favour of 
the advective-acoustic cycle as the origin of SASI. 

The difference at large shock radii between the 
advective-acoustic and the extracted timescale does not 
prevent the argument to be conclusive for several rea- 
sons. First, as compared to the large discrepancy with 
the purely acoustic time, this disagreement is quite slight. 
Second, it can easily be interpreted in the context of an 
advective-acoustic cycle by adjusting the effective cou- 
pling radius. Indeed the advective-acoustic time has been 
computed between the shock and the PNS surface, and as 
such should be considered as an upper limit. By contrast, 




(r sh -r.)/r. 



100 



Figure 5. Extraction of a radial propagation time from the 
unstable I = 1 SASI modes (Equation JJ), short dashed 
lines) and comparison with the advective-acoustic time (Equa- 
tion (TJ, full line) and the purely acoustic radial time (Equa- 
tion J2t . long dashed line). The times extracted from the eigen- 
spectrum are all much longer than the acoustic time (by a 
factor 3.5 to 7), which rules out a purely acoustic origin of 
these unstable modes. By contrast, they are comparable to 
the advective-acoustic time, thus confirming the advective- 
acoustic cycle as the origin of the instability. All timescales 
are in units of r^/v^. 



it is impossible to reconcile in the same way the acoustic 
time with the measured timescales, because the acous- 
tic time is too short. Finally we note that the instability 
mechanism at large shock radii was anyway not much de - 
bated, because the WKB analysis of lFoglizzo et al.l l|2007f ) 
can be conclusively applied to higher frequency harmon- 



4 COMPUTATION OF PURELY ACOUSTIC 
EIGENMODES 

In this section, a method to compute purely acoustic 
modes is described. In this method, one separates the ad- 
vected perturbations just below the shock, and subtract 
them from the total perturbations at the shock. Let us 
emphasise that a WKB description of the acoustic waves 
is not needed, since we do not need to distinguish acoustic 
waves propagating up and down. As a result, the method 
is applicable to lower frequency modes. We stress that 
we do not make any assumption regarding the transverse 
propagation of acoustic waves. The WKB method is ap- 
plicable when the wavelength of acoustic waves is shorter 
than the scale of the gradients: ur^/c^ S> 1, while the 
separation of advected perturbations only requires that 
the advection wavelength be shorter than the gradient 
scale: 



wr B h/v B h > 1- 



(10) 



For example for r s h = 2.5r», the first three / = 1 SASI 
modes have uur s h/v s h ~ 4.5 — 9 — 13. We may thus expect 
that the advected wave can be extracted reasonably well 
for the second and third modes, but only marginally for 
the first mode. 
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to i xr sh /v sh acoustic modes "WXh SASI modes 




Figure 6. Growth rate of the purely acoustic modes (left panel) compared with the growth rate of SASI (right panel). I = 1 
modes are shown in black and I = 2 modes in blue. For each spherical harmonics index, the most unstable mode is plotted (in the 
case of SASI, several harmonics are shown with thiner line for higher frequency modes). 



Another assumption implicitly made here is that, 
just below the shock, the entropy is simply advected, 
which is true if the flow is adiabatic. The flow just below 
the shock is indeed usually rather close to adiabaticity. 
The quasi-adiabaticity of the postshock flow can be quan- 
tified by the dimensionless parameter rdS/ dr(r s h) (where 
S is a dimensionless entropy defined in iFoglizzo et al.l 
which ranges from 7.10 4 for a shock radius at 
r s h = 5r„, to ~ 0.06 for r s h = 2.5r*. 

In order to be able to subtract the entropy and vor- 
ticity perturbations (described by the variables SS and 
Sq defined in Appendix [A"| . one needs to determine the 
values of the other variables Sf and Sh associated to this 
a dvected wave. This ca n be done as in the WKB ana lysis 
of lFoglizzo et al.1 (|2007r i; lYamasaki fc Foglizzol i|200Sl ). by 
requiring that these perturbations are advected: 



d_ 

Or 



(Sf) = 



Sf, 



™6h. 



(11) 

(12) 



Using the differential system recalled in Appendix A to 
remove the radial derivatives (neglecting the non adi- 
abatic terms in Equations (|A13|I - (|A14() due the quasi- 
adiabatic assumption), Equations (|11|) and (|12p can be 
expressed by the following linear system of two equations: 



M 2 6h adv 



l <5/adv 



Sh* 



= ss- 



l + (7-DM 2 ]f-§, (13) 
Sq 



(14) 



Note that Equation (|13[) (signifying that 8f is ad- 
vected) is equivalent to requiring that the pressure per- 
turbation vanishes (compare Equations (|13[) and ()A7|) ). 
Thus, when the advected wave is subtracted from the per- 
turbations below the shock, the pressure perturbation is 
unchanged. This system of equations can be solved to 
obtain the following expression of the advected pertur- 
bations as a function of the amplitude of the entropy 
and vorticity perturbations (described by the variables 



SS and Sq): 

Sfzdv = 
<5/ladv — 



1 - M 
I- H 2 M 2 
1 

1 - /i 2 M 2 
+ (1-M 2 )^ 



— SS - Sq 

7 



{^-(1-M 2 ) - 1)SS 



(15) 



(16) 



In order to compute purely acoustic modes, the ad- 
vected wave is then subtracted from the perturbations 
at the shock to define new boundary values <5/ acs h and 

Shvcsh'- 



S f&csh 
S /^acsh 



Sfsh — Sfadv, 
5h s h — (5/ladv, 



(17) 
(18) 



where Sf s h and 5h s h are determined by the boundary con- 
ditions at the shock (Equations (|A16[) and (|A16[l ). <5/ a dv 
and Sh&dv are computed using Equations (|15|l and (|16[) . 
and the values of SS and Sq at the shock (Equations (|A18f) 
and (|A19[) ). Using these modified boundary conditions, 
the acoustic modes can t hen be computed num erically as 
described for example in IFoglizzo et ail (|2007l ). 

Figure [U] compares the growth rate of the most un- 
stable acoustic modes (left panel) and SASI modes (right 
panel). The difference is striking: the purely acoustic 
modes are all stable for all shock radii in the explored 
range r s h = 2 — 100r» , while SASI modes have significant 
growth rates. This shows that the purely acoustic cycle 
cannot account for SASI and that the advected wave is 
essential to the instability. The most "unstable" acoustic 
mode (actually the least damped) for each I is the low- 
est frequency one. The discussion of the purely acoustic 
frequency spectrum is delayed to Section \5. 2. 2 1 

In order to quantify the effect of non-adiabaticity at 
the shock on the acoustic modes, we have used boundary 
conditions which either take into account or neglect the 
cooling processes in the computation of SS and Sq. The 
associated change in the frequency and growth rate of the 
modes is barely noticeable (less than 2% for any shock 
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radius), confirming the validity of our quasi-adiabatic ap- 
proximation. 

The error due to the high frequency approximation 
can in turn be quantified in the following way. Noting that 
there is some arbitrariness in the choice of the variables 
that should be advected, one can require another variable 
to be advected and compare the result with the previous 
calculation. For example requiring the perturbation of 
radial velocity Sv r /v r to be advected instead of Sh, we 
get a different definition of the advected perturbations: 



Sfudv 



6h B 



1 



1 + ei -M 2 (fi 2 + £2) 
-(1 + ei - M 2 {1 + c a )) 



(f + e 1 )(l-M 2 ) — 

7 

(19) 



5q 



T^—- 5S -—^ (20) 



where £i and £2 are two parameters vanishing in the limit 
cor s h/v s h S> 1, and defined by: 



02 



VjJ 

V 1 a 1 2 2 

~Or{C -V 



(21) 
(22) 



As expected, the two prescriptions for the advected per- 
turbations are equivalent in the limit of vanishing ei and 
£2- These parameters are very small for the second and 
third I — 1 acoustic modes: ei ~ £2 ~ 0.03 and 0.015 
respectively for r s h = 2.5r» (note that these two modes 
have a higher frequency than the second and third SASI 
modes as discussed in Section [5221) • The high frequency 
approximation is thus well justified and the two prescrip- 
tions are in reasonable agreement: the damping rates 
agree within 10 — 15% for the second mode, and < 5% 
for the third mode. These two modes are enough to con- 
clude that SASI is not a purely acoustic instability since 
the second and third I — 1 SASI modes are unstable. The 
approximation is slightly less good for the fundamental 
acoustic mode: in this case £1 ~ £ 2 ~ 0.1. The damping 
rates in the two prescriptions differ by ~ 30%, but the 
modes are still definitely stable whatever the prescrip- 
tion. 



5 DETAILED ANALYSIS OF THE 

ADVECTIVE-ACOUSTIC AND PURELY 
ACOUSTIC FREQUENCY SPECTRA 

In this section, we study in detail the frequency spectrum 
expected from advective-acoustic and purely acoustic cy- 
cles. In the sub section I5TT1 we study a simple toy model 
(Foglizzo 2009), which allows an analytical study of the 
frequency spectrum. In the following subsection l5.2.1l we 
then give an approximate generalisation to the more real- 
istic model of SASI used in Sections[3]and[4] This analysis 
confirms the validity of the argument presented in Sec- 
tionO Furthermore it allows a direct comparison between 
analytical expectations and the frequency spectrum of 
SASI. This direct comparison leads to the same conclu- 
sion about the instability mechanism, in favour of the 
advective-acoustic cycle. 



5.1 Planar toy model 

5.1.1 Presentation 



We consider the toy model introduced bv lFoglizzol l|2009h . 
It consists of a planar shock lying above a zone of de- 
celeration imposed by an external potential step. The 
potential step is localised within a length scale of H^, 
and is separated from the shock by a region of size H 
where the flow is uniform. The gas is assumed to be per- 
fect, adiabatic (with an adiabatic index 7 = 4/3), and 
non-self-gravitating. The direction perpendicular to the 
shock is called z, and referred to as vertical. Modes have a 
transverse structure along x (parallel to the shock), which 
is referred to as the horizontal direction. The domain is 
a box of horizontal size L. After a Fourier transform in 
time and horizontal direction x, the eigenmodes are char- 
acterised by their complex frequency uj and the number 
n x of horizontal wavelengths in the width of the box. The 
horizontal wave number is then: 



27m-, 



(23) 



The reader is referred to iFoglizzd (2009) for more details 
about this toy model, the equations governing it, and the 
numerical method used to solve them. In the remaining 
of this subsection, the numerical calculations have used 
the parameters: L — 6H, A4i = 00, H\r / H = 0.01, and 
cfJc 2 out = 0.75. 

Cycle constants Q and R describing the advective- 
acoustic and purely acoustic cycle can be defined as in 
Section [2] (but here without any approximation). An 
eigenmode with a complex eigenfrequency oj must then 
verify Equation ([5]), and contains both cycles in some 
proportion. In order to separate the two cycles and study 
their properties independently from each other, we follow 
IFoglizzd (2009) in considering advective-acoustic modes 
verifying: 



Q{u) = 1, 



(24) 



which corresponds to a situation where the acoustic wave 
propagating downward would be artificially removed. 
Similarly, purely acoustic modes follow: 



R{uj) = 1, 



(25) 



which signifies that the advected entropy- vorticity wave 
is ignored. This is equivalent to the method described in 
Section [4] to compute pur ely acoustic modes. 

As in IFoglizzd (2009). we define new cycle constants 
Qo and Ro that do not include the propagation of the 
waves between the shock and the potential step. These 
are defined as: 



Q 
R 



Qoe 



i(fcJ dv -fcJ)H 



R Q e^- k * )H , 



(26) 
(27) 



where fcj dv and kf are the vertical numbers of the ad- 
vected wave, and of the acoustic waves propagating down 
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(+) and up (- 



kf = - 



Vz 



= 1 - 



cl-M 2 

k 2 x c 2 (l~M 2 ) 



(28) 
(29) 
(30) 



5.1.2 Eigenfrequencies of the advective- acoustic cycle 



Equations (|24[1 and (|26[l are combined to obtain the con- 
dition verified by an advective-acoustic mode: 



1+uM 

Q oe — t+aT = 1, (31) 

where r aac is the duration of a one dimensional advective- 
acoustic cycle (k x = 0, = 1): 

^ E. — — 

r aac = — + c( ^_ M ~ ) - \ v \(i_ M y 

Equation [31] can be solved approximately in the regime 
where the growth rate is low (u>i -C tJr) and where the 
acoustic wave is propagating (/i 2 > 0). In this case, 
the modulus of this comple x relation gives access to the 
growth rate (|Foglizzoll2009l 'l: 



(32) 



l + M log \Q\ 

7 fJ>r • 



(33) 

fir + M T aac 

In a similar way, the frequency may be obtained from the 
argument of the same relation, which gives an informa- 
tion on the wave phase: 



$(Qo) + Re OJT; 



1 + 



l + M J 



2n z 7r, 



(34) 



where $(Qo) is the phase of Qo, and describes the phase 
shift taking place during the two coupling processes. 
Equation (|34p gives a second order polynomial in u) r , 
which physical solution is: 



1 



1-M 

-M [(n 



(n z - $/27r)Lj a 



-<&/2tt) 



- nW^il - M) 2 } l7 (3f ) 



where we defined two frequencies typical of the advective- 
acoustic cycle and of the horizontal propagation: 

2tt 



"7~aac 
kx C. 



(36) 
(37) 



Note that in principle $ depends on the frequency, there- 
fore Equation (|35[) gives explicitly the frequency only if 
we know the value of $. In the following of this subsec- 
tion, we suppose $ = which turns out to be a good ap- 
proximation in this toy model. This is not always true in 
the spherical model as will be discussed in Section [5.2.21 
Figure [7] compares this analytical result with eigen- 
modes computed numerically. The agreement is very 
good both with advective-acoustic modes and with the 
full modes (i.e. the modes that verify Equation ©), 
in the regime of propagative acoustic waves. The left 
panel shows that the approximations made are valid, and 
that there is no significant phase shift. The right panel 



shows that the frequency of the full modes is dictated by 
the advective-acoustic c ycle, as could be expected from 
the results of iFoglizzol (2009). One can nonetheless no- 
tice that the purely acoustic cycle has an influence on 
the growth rate of the modes (esp ecially those close to 
marginal stability, (|Foglizzo|[2009l )'). since the unstable 
modes are not exactly the same on the two panels. The 
origin of the stable n x = 1 mode at very low frequency, 
which does not coincide with any curve, is rather uncer- 
tain (it lies in the regime of evanescent acoustic waves, 
which is not described by Equation (I35|) ). 

From Figure [7] and Equation (|35|) , one may notice 
that two successive modes sharing the same transverse 
structure are separated by a frequency difference that 
is close to the frequency of a vertical advective-acoustic 
cycle cj aac . This is exact in the case of one-dimensional 
modes, but only approximative for n x 7^ 0. Indeed a 
transverse advective-acoustic mode has a slightly higher 
frequency than the corresponding ID mode (by a factor 
at most 1/(1 — M) for a propagative mode). This higher 
frequency may seem counterintuitive at first sight, be- 
cause the geometrical path along which the waves prop- 
agate during one cycle is longer. This longer path is the 
reason why the duration of the cycle increases with incli- 
nation in Equation (|33|l giving the growth rate. However 
this duration corresponds to the group velocity of acous- 
tic waves, while the frequency is determined by a phase 
relation. The cycle time relevant to the frequency is thus 
given by the phase speed. The vertical phase speed of 
an acoustic wave increases with inclination, because the 
wave fronts get closer to the vertical direction. This ex- 
plains the higher frequency of non radial modes. 



5.1.3 Eigenfrequencies of the purely acoustic cycle 

The frequency of the purely acoustic modes can be ob- 
tained by a very similar method to that used for the 
advective-acoustic modes. For this purpose, one combines 
Equations (l2"5"Tl and P7)l to obtain: 



7-) u 

R e 



= 1, 



(38) 



where r acz is the duration of a one-dimensional purely 
acoustic cycle (k x = 0): 

H H 2H 



+ 



c{l + M) c(l-M) c(l-M 2 )' 



(39) 



By taking the argument of this complex equation, we 
obtain: 



$(i?o) + 7?e(r acz ^ 2 -n 2 Lu 2 cx (l - M 2 )) = 2n z 7r. (40) 

from which one can deduce the following expression for 
the frequency of purely acoustic modes: 



Lu r = [{n z - $/2vr) 2 tj acz + n 2 x ujl cx (l - M 2 )] 



1/2 



(41) 



where aj acz is the frequency associated to a vertical acous- 
tic propagation: 

Lu acz = . (42) 

T a c2 

Equation (|4ip is analogous to the frequency of the acous- 
tic modes of a rectangular box with periodic boundary 
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Figure 7. Eigenspectrum of the advective-acoustic cycle in the planar toy model of iFoglizzol (|200gh : oscillation frequency as 
a function of n x (the number of transverse wavelength in the width of the box). Symbols represent the numerically computed 
eigenmodes for the advective-acoustic cycle only (left panel) and the full toy model (right panel). The colour and shape of the 
symbols reveal minor differences in their stability properties: turquoise diamonds for stable, purple circles for unstable, and red 
squares for the most unstable mode for a given value of n x . Continuous lines represent the analytical solution given by Equation l|35ll 
and assuming 4> = (the curves are restricted to the regime of propagating acoustic waves, where the analytical expression applies). 
The analytical expression is in very good agreement with the numerical modes in both panel, as would be expected from the fact 
that the advective-acoustic cycle is driving the instability in this toy model. 



co/co purely acoustic modes co/co full modes 

acz ' acz 




Figure 8. Same as Figure[7]but for the purely acoustic cycle. The analytical expression (Equation 1411 1 . with <E> = 0) agrees with 
the purely acoustic modes (left panel), but not with the full modes (right panel) confirming that these are not dominated by the 
purely acoustic cycle. 



conditions in x and reflexive boundary conditions in z (it 
is only slightly modified by advection): 



2ttc 



(2Hf 



(43) 



Thus as argued in Section[3] if the number of wavelengths 
in the transverse direction (n x ) increases, the frequency 
increases typically by cj ac:c corresponding to a propaga- 
tion time in the transverse direction. Similarly if the num- 
ber of wavelength in the vertical direction (n z ) increases, 
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the frequency increases typically by uj RC z corresponding 
to an acoustic propagation along the vertical direction. 

Equation (|41[) (with $ = 0) is compared with the 
numerically computed eigenmodes in Figure [8] It agrees 
with the purely acoustic modes (all stable), showing here 
again that the approximation of small growth rate holds 
and that the phase shift is not very significant. By con- 
trast, the frequency of the full modes do not match the 
purely acoustic prediction. This confirms that the unsta- 
ble modes are not dominated by the purely acoustic cycle, 
and that the study of the frequency spectrum enables to 
discriminate between the two mechanisms. 

Although the purely acoustic cycle does not deter- 
mine the mode frequencies, it can be observed from Fig- 
ure [8] that the acoustic timescale influences the modes 
stability. Indeed, the full modes in the vicinity of an 
acoustic branch are more likely to be unstable than the 
other modes. Furthermore, the most unstable modes (for 
a given n x ) are all close to the branch n z = 0. These ob- 
servations can be interpreted by two phenomena. First, 
a constructive interference between the purely acoustic 
and advective-acoustic cycles takes place when the mode 
frequency (determined by the advective-a coustic cycle) i s 
close to that of a purely acoustic cycle (Foglizzo 2002). 
This destabilising effect is more pronounc ed for modes 
close to marginal stability ijFoglizzd 120091 ). Second, the 
growth rate associated with the advective-acoustic cycle 
alone is maximu m when the aco ustic propagation is close 
to horizontality (Foglizzo 2009), which happens at the 
frequency of the fundamental acoustic branch n z = 0. 



nar model oj aac , o; ac , , and uj. 



5.2 Spherical model 

5.2.1 Approximate generalisation of the analytical 
results 

The description of waves propagation given in the previ- 
ous subsection is exact in a uniform flow, and could be 
generalised to a non uniform flow in spherical geometry 
with the use of the WKB approximation. In this frame- 
work, the terms of the form e lk * H (in Equations (126 p 
and I27[) should be replaced by terms of the following 



form: 



f k r dr 



However, as our goal is to discriminate 



the mechanism at work in the regime where the WKB 
approximation does not hold, we will not develop rigor- 
ously such a description. We simply give an approximate 
generalisation, from which one may expect the right or- 
der of magnitude but not an exact value. This is however 
sufficient to distinguish the two mechanisms as noted in 
Section [3] 

The timescales associated to a radial advective- 
acoustic cycle, a radial purely acoustic cycle, and a trans- 
verse acoustic propagation have been defined for a spher- 
ical model in Section 12.31 (Equations (UJ-©). One can 
define from these timescales, typical frequencies which 
are a direct generalisation of the frequencies of the pla- 



aac , ^ac2 , 

2tt 

^aac — j 
7~aac 

2tt 

^"acr — j 
Tacr 

f . 2-k c 
Wa,c<j>{r) = = -. 

7~ac(f> i 



(44) 
(45) 
(46) 



We propose to estimate the frequency spectra of 
advective-acoustic and purely acoustic cycles by directly 
introducing these typical frequencies in Equations (|35[) 
and (|41l) (for this purpose n^w| should be replaced by 



acc/> / 



. This procedure suffers from different sources 



of uncertainty: the formula itself is approximate, the 
phase shift is a priori unknown, and the typical frequen- 
cies c^aac, i^acr, and u) ac <p are uncertain respectively due 
the lack of knowledge of the coupling radius, the turn- 
ing point of acoustic waves, and the radius at which it 
should be evaluated. When comparing to SASI modes in 
Section T5. 2. 21 we will use a range of values for the radius 
where a; ac and M are evaluated in order to illustrate 
this uncertainty. 

The above formulation can be used to assess the va- 
lidity of the argument developed in Section [3] An inspec- 
tion of Figures [7] [H] [5] and [TO] confirms visually that the 
difference between two successive modes that share the 
same transverse structure corresponds approximately to 
a radial propagation frequency: oj aac for the advective- 
acoustic cycle, and c^acr for the purely acoustic cycle. 
This visual observation can be recovered from Equa- 
tion (|35p for the advective-acoustic cycle, if one assumes 
1(1 + l)^ac0 < n?Wa ac : 



Aw 



I - 



2tt 



M)l(l + lVac0 



2(n r - $/27r)(n r + I - $'/2tt)w| 



+ 



(47) 



where we have defined Aco = uj r (n r + 1) — u) r (n r ), 
$ = $(n r ), and $' = $(n r + 1). Similarly for the purely 
acoustic cycle, Equation (|4ip with the assumption that 
1(1 + l)cu1 c( p <C n^Wacr leads to the following relation: 



Aw 



2tt 
(1- 



■M 2 )l(l + l)u&+ 



2(n r - $/27r)(n r + I - $'/2ir)w2 



+ 



(48) 



Thus we expect the timescale extracted from the eigen- 
frequencies in Section [3] to match one of the radial times 
(depending on the instability mechanism), for sufficiently 
large values of n r and if the two modes have the same 
phase shift. The deviation at low n r can be estimated 
from the second term in the right hand side of Equa- 
tions (|47[) an d (|48[) . For example for the parameters 
r s h = 2.5r», I — 1, n r = 1, it is already less than ~ 10% 
for both the advective-acoustic and the purely acoustic 
modes (note that for the purely acoustic modes n r — 1 is 
the second mode) . This confirms the validity of the argu- 
ment given in Section |3j and is roughly consistent with 
the fact that the three extracted times differ by at most 
15 — 20%. We note that a difference in the phase shift 
between two successive modes could pollute the measure 
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i r sh h purely acoustic modes 

50 I ' 1 




Figure 9. Comparison of the frequency spectrum of the 
purely acoustic eigenmodes (computed as described in Sec- 
tion [4t with the analytica l estim ate, in the spherical model 
of iBlondin fc Mezzacappal 1200$) . The frequency in units of 
^sh/ r sh is represented as a function of the spherical harmonics 
index I. Symbols represent eigenmodes with the same colour 
and shape convention as the previous figures. Full lines corre- 
spond to analytical estimates in which Jv[ and w ac are cal- 
culated at the radius r = (r^ — r*)/2, while dashed lines use 
r = r B h and r = l.lr*. In order to obtain a decent agreement, 
the phase shift of the harmonics (n r > 0) had to be adjusted 
to the value <J>/(2tt) = 0.4. 

of the radial time, but it is not expected to consistently 
change this measure for different mode duplets. 



Figure [10] compares the frequency of SASI eigen- 
modes with the predictions from the advective-acoustic 
cycle (left panel) and the purely acoustic cycle (right 
panel). The advective-acoustic cycle prediction agrees 
with the frequency of all the non radial modes, with a 
precision of ten to a few tens of percents. On the other 
hand, the purely acoustic cycle approximately predicts 
the frequency of the most unstable modes for each I, but 
significantly overestimates that of other modes. This is 
consistent with the result of Section [3] that the acoustic 
radial time is significantly shorter than that extracted 
from the eigenspectrum. It is also consistent with the fact 
that the SASI eigenspectrum is much more densely popu- 
lated than the acoustic eigenspectrum (compare Figure [9] 
and Figure 1101 remembering that the frequency scale is 
different). 

Similarly to the toy model, the most unstable SASI 
modes for each spherical harmonics index I are close 
to the purely acoustic branch n r — 0. This can again 
be interpreted by a higher efficiency of the advective- 
acoustic cycle, as well as a constructive interference with 
the purely acousti c mode. It is also probabl y linked with 
the observation of IFernandez fc Thompson! (|2009bh that 
the growth rate of a mode is maximum when the advec- 
tion time coincides with the azimuthal acoustic time. 

Finally, we observe that the (stable) radial modes 
are systematically shifted with respect the theoretical 
curve. This can be explained by a phase shift taking 
place during the advect i ve-acou stic cycle, as shown by 
IFernandez fc Thompson! l|2009al ). Indeed, a phase shift of 
$ = n then gives the following frequency spectrum for 
the I = modes: 

Lo r = (n r - l/2)cj aac , (49) 

which is in much better agreement with the numerical 
eigenmodes. 



5.2.2 Comparison with SASI and purely acoustic 
eigenfrequencies 

Figure [9] compares the frequency spectrum of purely 
acoustic modes with the analytical estimate. The method 
described in Section [4] gives two kinds of modes: the 
acoustic modes in which we are interested, as well as 
modes due to a secondary advective-acoustic cycle with 
the entropy being created by cooling processes acting 
on the acoustic wave. This last category of modes are 
clearly distinguished and significantly more stable than 
the purely acoustic modes, and therefore we do not show 
them in this discussion of the purely acoustic modes. Con- 
trary to the case of the planar toy model, it is necessary to 
invoke a phase shift so that the analytical estimate agrees 
with the numerically obtained eigenfrequencies. We find 
that a phase shift of $/(2tt) = 0.4 for all modes with 
n r > 0, and $ = for all modes with n r = gives a 
reasonable match to the eigenspectrum (Figure [9]). The 
origin of this phase shift remains obscure. Nevertheless, 
the comparison confirms that the analytical estimate ac- 
counts for the frequency of acoustic modes within the 
uncertainties discussed in the previous subsection and il- 
lustrated by the dashed lines in Figure [9] 



6 CONCLUSION 

Two new arguments have been presented which allow us 
to conclude that SASI cannot be explained by a purely 
acoustic mechanism, and should thus be interpreted as 
an advective-acoustic cycle. Contrary to previous works, 
these two arguments do not rely on a WKB analysis of 
the acoustic waves, and are therefore valid for realistic 
values of the shock radius. 

First we proposed a method to extract a radial prop- 
agation timescale from the frequency difference between 
two successive SASI modes. This timescale agrees fairly 
well with the advective-acoustic time, but strongly dis- 
agrees with the radial acoustic time. This discards the 
purely acoustic cycle as the driving mechanism of SASI, 
and argues in favour of the advective-acoustic cycle. A 
detailed analysis of the frequency spectrum provided in 
Section [S] confirmed the validity of this argument. 

Second, we described a new method to compute 
purely acoustic modes by artificially subtracting the ad- 
vected wave created by the shock. All these acoustic 
modes were found to be stable, thus showing that the 
advected entropy-vorticity wave created by the oscilla- 
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advective-acoustic predictions co r r s/ V sh purely acoustic predictions 




Figure 10. Comparison of the frequency of SASI eigenmodes with the estimates for an advective-acoustic cycle (left panel) and 
a purely acoustic cycle (right panel). The frequency in units of v s ^/r s ^ is represented as a function of the spherical harmonics 
index I. Symbols represent eigenmodes with the same colour and shape convention as the previous figures. Full lines correspond 
to analytical estimates in which M and w ac are calculated at the radius r = (r a ^ — r»)/2, while dashed lines use r = r s ^ 
and r = l.lr*. As previously, the advective-acoustic frequency estimate is represented only in the regime of propagating acoustic 
waves (at the radius r used for the estimate). The prediction from the advective-acoustic cycle is roughly in agreement with all 
the non radial modes. The purely acoustic cycle approximately predicts the frequency of the most unstable modes (for a given I), 
but strongly overestimates the frequency of higher frequency harmonics. 



tion of the shock plays a central role in the instability 
mechanism. 

As a side product, we provided an analytical descrip- 
tion of the frequency spectra of the advective-acoustic 
and purely acoustic cycles . This descripti on was checked 
in the planar toy model of iFoglizzol (|2009f ) and the model 
of lBlondin &: Mezzacappal (|200rj ). The similarity between 
the frequency spectra in these two models, supports the 
relevance of the toy model for core collapse supernovae. 
Finally, we note that unexplained phase shifts during the 
coupling processes were observed in the second model. 
This calls for a better understanding of the coupling pro- 
cess in a cooling layer. 



APPENDIX A: DIFFERENTIAL SYSTEM 
AND BOUNDARY CONDITIONS WITH THE 
NEW VARIABLES OF YAMASAKI & 
FOGLIZZO 



Following lYamasaki fc Foglizzol (|2008b . we adopt the 
variables SS, Sq, Sf, Sh defined by: 



ss 



2 Sc Sp 
P 

dr 





7 — 1 c 


Sq = 


\J PVr 


Sf = 


V r Sv r H 

7 - 


Sh = 


5v r Sp 
+— , 



-cSc — Sq, 



(Al) 
(A2) 
(A3) 
(A4) 
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where S denotes the Eul erian perturba t ion. A s compared 
to the variables used bv lFoglizzo et al. I l|2007l ). SS and Sh 
are unchanged, the new variable Sq is related to SK by 
SK = l(l + l)5q, and Sf is slightly changed to incorporate 
the perturbed heating Sq. The usual physical quantities 
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can be expressed as a function of these new variables as: 



Sv r 1 

~ ~ 1-M 2 

8p 1 

~p 1-M 2 

8P 1 

7P = 1 - M 2 

Sf + Sq 



5h + 5S - s 4- s 4 



(A5) 



-M 2 5h-SS + 



Sf + 5q 



(A6) 



M 2 Sh-{l + ( 1 -l)M 2 ) — 



+ 



Sc 2 



7-1 fSf + 5q 



1-M 2 



(A7) 



M 2 Sh - M 2 8S . (A8) 



The velocity in non radial directions has a different spa- 
tial structure and is expressed as a function of angular 
derivative of the spherical harmonics in the following way: 



Svg — 



Svd, = 



Sf BYT 
iujr 89 ' 

Sf ovr 



(A9) 
(A10) 



The advantage of this new formulation is that the 
differential system describing the perturbations becomes 
particularl y compact (note tha t it is strictly equivalent 
to that of iFoglizzo et all (120071 ) ): 



dSf 
dr 



dSh 
dr 



dSS 
dr 
dSq 
dr 



IU)C 



1 M 2 Sh - M 2 ^ 
v r (l - M 2 ) \ c 2 

+ [l + (7-l)>! 2 ]--^}, (AH) 
7 c J 



V r (l -M 2 ) \ C 



11 Sf-M 2 Sh 



-8S- 



Sq 



* w ss + s ( ^ 

v r V Pv, 



— Sq + S[ , 

Vr \P v r / 



where is defined as in IFoglizzo et ail (|2007l ) by: 



l(l + l)c 



-(1-M 2 ), 



(A12) 
(A13) 
(A14) 

(A15) 



and the exp licit expression s of 5 (£/pv r ) and 8(C/Pv r ) 
are given in IFoglizzo et al.l (|2007h (Equations (A7) and 
(AS))- 

The boundary conditions at the shock are obtained 
by imposing the Rankine-Hugoniot relations for the per- 
turbed quantities, which are written as: 



5/ s h = iujVr,lA( I 1 



"r,sh 



Sh sh = AC ( 1 

^r,sh 

SSsh .IjJVrA \ j- I 1 

— = *— 2— A C 1- 



V r ,l 
V r ,l 



(A16) 
(A17) 



Sq B 



\ Vr,l J Psh«r,sh C sh 

v r, s h \ AC ( 2u r> i« P , S h d$ 

«r,l / c; 
£sh — £\ 

Psh^r,sh 



C 2 



-AC, 



^J.(A18) 
(A19) 



where the subscripts 'sh' and T' refer to the values just 
below and above the shock, respectively; AC is the radial 
displacement of the shock surface. 
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